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Abstract 


The  niunerical  solutions  of  two  related  variations  of  a  sequential  version  of  a  form  of 
the  bioequivalence  problem  was  presented  in  a  report  by  Hwang  (1991).  In  that  report 
he  referred  to  our  unpublished  results  bounding  the  solutions  and  providing  asymptotic 
expansions.  The  present  report  has  two  major  functions.  One  is  to  derive  and  amplify  these 
results,  and  incidentally  to  correct  an  error.  The  second  is  to  gather  together  in  one  place, 
and  with  relatively  little  of  the  abbreviation  characteristic  of  previous  publications,  many 
of  the  details  that  are  useful  in  deriving  the  asymptotic  expansions,  e.g.,  see  Breakwell  and 
ChemofF  (1964),  Chemoff  (1965a,  1972),  and  Chemoif  and  Petkau  (1981). 

In  Hwang’s  report  he  presented  the  Bayesian  decision  theoretic  approach  in  which 
the  problem  is  related  to  a  stopping  problem  involving  Brownian  motion,  a  stopping  cost 
represented  by  a  single  cost  parameter  c,  and  an  initial  point  (y,  s)  depending  on  the 
prior  normal  distribution  of  an  unknown  parameter  and  other  known  parameters  of  the 
problem.  The  solution  of  the  problem  is  represented  by  dividing  up  the  set  {(y,  s)  :  s  >  0} 
into  a  continuation  set  C  and  a  stopping  set  S. 

AMS  1991  subject  classification.  Primary  62L10,  secondary  62P10. 

Key  words  and  phrases.  Bioequivalence,  sequential  anlaysis,  Baye&  risk.  Brownian  motion, 
free  boundary  problems,  optimal  stopping,  asymptotic  expansions. 


1.  Introduction 


The  n\imerical  solutions  of  two  related  variations  of  a  sequential  verrion  of  a  form  of 
the  bioequivalence  problem  was  presented  in  a  report  by  Hwang  (1991).  In  that  report 
he  referred  to  our  impublished  results  bounding  the  solutions  and  providing  asymptotic 
expansions.  The  present  report  has  two  major  functions.  One  is  to  derive  and  amplify  these 
results,  and  incidentally  to  correct  an  error.  The  second  is  to  gather  together  in  one  place, 
and  with  relatively  little  of  the  abbreviation  characteristic  of  previous  publications,  many 
of  the  details  that  are  useful  in  deriving  the  asymptotic  expansions,  e.g.,  see  Breakwell  and 
Chemoff  (1964),  ChemofF  (1965a,  1972),  and  ChemolF  and  Petkau  (1981). 

In  Hwang’s  report  he  presented  the  Bayesian  decision  theoretic  approach  in  which 
the  problem  is  related  to  a  stopping  problems  involving  Brownian  motion,  a  stopping  cost 
represented  by  a  single  cost  parameter  c,  and  an  initial  point  (y,  s)  depending  on  the 
prior  normal  distribution  of  an  unknown  parameter  and  other  known  parameters  of  the 
problem.  The  solution  of  the  problem  is  represented  by  dividing  up  the  set  {(y,  s)  :  s  >  0} 
into  a  continuation  set  C  and  a  stopping  set  S. 

In  Section  2,  we  shall  describe  the  formulations  of  the  bioequivalence  problems,  the 
corresponding  stopping  problems  and  the  relevance  of  the  free  boundary  problem.  In 
Section  3  we  describe  the  results  of  Hwang  and  how  our  theoretical  results  related  to 
these.  In  Section  4,  we  cumulate  details  about  separable  solutions  of  the  heat  equation. 
These  will  serve  as  a  repository  for  future  reference  and  many  of  these  will  be  relatively 
unnecessary  for  this  particular  problem. 

In  Section  5  we  will  derive  the  bounds  and  in  Section  6,  the  asymptotic  expansions  will 
be  calculated  formally.  No  attempt  will  be  made  to  repeat  the  type  of  proof  of  Chemoff 
(1965a)  that  these  expansions  are  true  asymptotic  expansions. 

For  these  particular  problems,  Hwang’s  computer  program,  modeled  on  Chemoff  and 
Petkau  (1984)  provides  a  very  good  numerical  solution,  and  the  bounds  and  expansions  are 
of  little  more  than  academic  interest.  Nevertheless  these  theoretical  results  were  derived 
first,  and  provided  a  qualitative  view  of  the  solution  which  seemed  quite  strange  at  the 
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time.  Without  that  view,  it  may  have  been  difficult  to  produce  a  good  computation.  In 
one  particular  case,  Hwang  shows  when  the  accuracy  of  the  numerical  procedure  has  to  be 
improved,  as  it  readily  can  with  more  computer  time,  in  order  to  match  better  with  some 
of  the  asymptotic  expansions. 

Most  of  the  boimds  in  Section  5  apply  relatively  simple  principles  compared  to  those 
pioneered  by  Bather  (1962,  1970)  and  illustrated  in  ChernofF  (1972).  In  spite  of  the  length 
of  this  msmuscript,  the  amount  of  time  the  authors  spent  in  their  collaboration  leading  to 
these  theoretical  results  was  relatively  small,  because  it  was  felt  that  enough  was  known 
to  facilitate  the  numericad  calculation.  However  a  miscEdculation  at  the  time  failed  to 
indicate  the  presence  of  another  term  of  theoretical  interest,  which  is  still  unaccounted  for 
and  which  will  be  discussed  in  Section  6. 

2.  Problem  Formulation  and  Summary  of  Results 

There  are  several  alternative  formulations  of  the  bioequivalence  problem.  Some  time 
ago  Bather,  ChernofF  and  Petkau  (1988)  presented  the  following  formulations.  In  his 
description  of  the  formulation,  Hwang  referred  to  work  by  Metzler  (1974),  O’Quigley  and 
Baudoin  (1988)  and  Selwyn  et.  al.  (1981). 

Let  X  be  a  normally  distributed  variable,  with  unknown  mean  /i  and  known 
variance  <t^,  representing  the  difference  between  the  outcomes  of  two  treatments  on  a 
matched  padr  of  subjects.  One  of  these  treatments  is  the  staindard  aind  the  other  is  a  new 
treatment  which  is  hoped  to  have  matching  pharmacological  properties.  As  observations 
on  X  cumulate,  the  decision  is  made  to  continue  or  stop  sampling.  When  sampling  is 
stopped,  the  decision  is  made  to  declare  difference  or  equivalence.  Presumably  if  difference 
is  declared,  the  search  for  a  new  equivalent  formulation  is  repeated. 

The  costs  involved  are  that  of  having  to  start  over  if  nonequivalence  is  announced, 
that  of  having  to  accept  the  mean  /x,  which  may  be  different  from  0,  when  equivalence 
is  announced,  and  a  cost  of  sampling.  The  cost  of  starting  over  is  some  constant  ko 
independent  of  fi  since  fi  is  simply  a  characteristic  of  the  candidate  drug  which  is 
discarded.  The  cost  of  accepting  bioequivalence  should  depend  on  how  large  is.  We 
will  consider  two  versions.  In  the  first  this  cost  is  kifi^  and  in  the  second  it  is  fc2|;x|. 
Finally,  we  shall  assume  that  the  cost  of  sampling  is  lineair  in  the  sample  size  n,  i.e.,  c^n 
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for  n  pairs  observed.  Implicit  in  this  sequential  version  of  the  problem  is  the  assumption 
of  immediate  observed  response  to  treatment. 

Clearly  the  operating  chau’acteristics  of  any  procedure  applied  in  this  problem  will 
depend  on  the  unknown  mean  y..  The  concept  of  the  optimal  procedure  requires  some 
criterion  of  optimality.  By  accepting  a  Bayesian  formulation  with  a  prior  normal  distribu¬ 
tion  with  mean  and  variance  ,  we  establish  a  clear  cut  optimization  problem,  which 
can  in  principle  be  solved  directly  by  the  backward  induction  of  dynamic  programming. 

The  theoretical  study  of  this  problem  can  be  facilitated  by  converting  the  problem 
to  a  related  one  where  the  sum  of  the  first  n  i.i.d.  observations  Xi,  A’2, . . . ,  JCn  is 
replaced  by  Brownian  Motion  X(t)  with  unknown  drift  y  and  known  variance  per 
unit  time.  The  original  discrete  time  problem  may  be  regarded  as  a  special  version  of  the 
continuous  time  problem  where  stopping  is  permitted  only  at  certain  discrete  times.  One 
outcome  of  this  continuous  time  formulation  is  that  the  class  of  problems  can  be  reduced  to 
a  one  parameter  family  of  problems,  the  numerical  solution  of  which  can  be  approximated 
arbitrarily  well  by  techniques  outlined  in  Chernoff  and  Petkau  (1986).  Another  outcome  is 
the  relation  of  the  optimization  problem  to  the  solution  of  a  free  boundary  problem  (FBP) 
involving  the  heat  equation. 

In  the  discrete  time  Bayesian  version,  the  posterior  probability  distribution  of  y  after 
n  observations  is  given  by 

(2.1) 


y  _  yo<^o  ^  +  (-^1  Xn  ^ 

-I-  71(7"^ 


(2.3) 

In  the  continuous  time  version 

(2.4)  C{y\X{t'), 0  <  <  t)  =  N{Y{s),s) 

where 

(2.5)  Y{s)  =  {yoGo^  X{i)G-^)/{G^^  +  tG-^) 


=  {CTo  ^  +n(7  ^)“^  . 


where 

(2.2) 

and 


Ciy\X„X2,....X„)  =  NiYn.s,,) 
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and 


(2.6)  s  =  (<T^^  +ta  ^  . 

Moreover,  it  can  be  shown  that  y(5)  is  a  Brownian  motion  with  zero  drift  in  the  —s 
scale  originating  at  (t/c^o)  = 

(2.7)  £^[dy(5)]  =  0 

(2.8)  Var(dy(s)]  =  -ds 

Note  that  as  t  increases  from  0  to  oo,  s  decreases  from  sq  to  0. 

The  posterior  risk  in  stopping  at  y(s)  =  y,  is  given  in  the  first  version,  which  we 
shall  refer  to  as  Problem  1,  by 

min(*:o,-E[/^^|y(^)  =  yj)  =  min(fco, /:i(y^  +  «))• 

In  the  second  version  or  Problem  2,  the  posterior  risk  can  be  czdculated  to  be 

m-(fco,JE:[lpliy(s)  =  y])  =  min(4, 


where 

(2.9)  Gu{a)  =  2{<i>ia)  +  a[$(o)  -  1/2]} 

and  <i>  and  $  are  the  steindard  normal  density  and  cumulative.  Thus  the  stopping  risk 
in  Problem  1  is  given  by 

du(y,s)  =  cot  +  min(Aro,  fci(y^  +  s)) 

(2.10)  =  co<T^s~^  +  min(A:o,  fci(y^  +  «))  - 

and  in  Problem  2  by 

(2.11)  d2i(y,  s)  =  +  min(A:o,  k2S^^'^Gu{ys~^^^))  ~  cocr^crj"^. 

For  a  stopping  cost  d{y,s)  there  is  a  corresponding  optimal  Bayes  Risk  p(y,  s),  given 
that  we  have  reached  y(5)  =  y,  and  proceed  optimally  thereafter.  But  then  it  pays  to 
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stop  if  p{y,s)  =  d{y,s)  and  to  continue  if  p(y,s)  <  d(y,s).  But  this  means  that  the 
optimal  stopping  time  5  corresponds  to  a  procedure  which  stops  when  (K(s),  s)  reaches 
a  stopping  set  S  =  {(y,s)  :  p(y,s)  =  That  is,  the  optimal  procedure  is  one  of  a 

class  of  procedures  which  is  defined  by  a  set  S  and  its  complement  C  =  S’^. 

Any  stopping  time  5,  optimal  or  suboptimal,  which  is  defined  by  a  subset  6'  of 
{(y,  s)  ;  s  >  0},  has  a  corresponding  Bayes  risk  given  by 

(2.12)  6(v,‘>)  =  £Ry(S),S)ll'(»)  =  if) 

which  satisfies  the  ordinary  boundary  problem 

(2.13)  ~  (yi"*)  ^  ^ 

and 

(2.14) .  6(y,  s)  =  d(y,  s)  for  (y,  s)  e  S  . 

The  extra  condition  that  determines  the  optimality  of  S  is 

(2.15)  by(y,s)  =  dy{y,3)  for  (y,s)  6  dS 

where  dS  is  the  boundary  of  S.  The  pair  {S,b)  which  satisfy  (2.12)  -  (2.14)  represents 
the  solution  of  the  free  boundary  problem  (FBP)  and  relatively  mild  sufficiency  conditions 
under  which  the  solution  of  the  FBP  represents  the  solution  of  the  optimum  stopping 
problem,  are  discussed  in  Chernoff  (1972).  To  paraphrase  those  conditions,  they  state  that 
if  the  stopping  rule  determined  by  the  solution  of  the  FBP  can’t  be  improved  trivially,  the 
stopping  rule  is  optimeil. 

Since  the  stopping  set  S  does  not  depend  on  the  initizd  value  of  (y,  s),  i.e.,  (po,  o'o)? 
we  see  that  the  set  S  represents  a  solution  which  applies  simultaneously  for  all 
The  remaining  parEuneters  cq,  ko,  ki  in  Problem  1  can  be  incorporated  into  one  essential 
parameter  by  use  of  a  transformation  of  the  form 

(2.15)  F*  =  aY,  s*  =  a^s 

which  leaves  fJ[dF*(s*)]  =  0  and  Var[dF*(s*)]  =  — ds*,  and  by  noting  that  the  solution 
S  of  the  FBP  is  not  altered  when  d  is  multiplied  by  a  constant  or  modified  by  the 
addition  of  a  solution  of  the  heat  equation. 
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Now  let  fci/a^fco  =  1,  i.e., 

(2.17)  a  =  (ki/koY^^  . 

Then 

(2.18)  duiy,s)  =  kodiiiy\s*)  -  cocr^aj'^ 
where 

(2.19)  di2{y,s)  =  cis~^  +  inin(l,y^  +  s) 
and 

(2.20)  Cl  =  co<T^ki/kl 

and  the  solution  of  the  FBP  or  optimization  problem  for  di2  can  be  converted  simply 
to  that  for  dn.  Note  that  we  now  have  a  FBP  with  only  one  parameter  ci  to  replace 
/io,a-2,co,fco  and  ki. 

A  similar  transformation  works  for  Problem  2.  Here  we  have,  with  a  =  k2/ko 

(2.21)  d2iiy,s)  =  kod22{y*,s*)  -  co^^crj'^ 
where 

(2.22)  d22iy,s)  =  C25“‘  +min(l,s^/^Gie(ys“^^^)) 
and 

(2.23)  C2  =  coc^klfkl 

Later  we  shall  see  computational  advantages  for  deriving  zisymptotic  expansions  by 
subtracting  special  solutions  of  the  heat  equation  from  di2  and  <£22- 

3.  Character  of  the  Solutions 

Because  of  the  symmetry  of  the  loss  functions  eis  functions  of  /x,  5,  the  optimal  S 
and  S‘^=C  are  symmetric  in  y.  These  sets  are  adequately  described  by  two  curves,  an 
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inner  and  an  outer  curve  given  by  yi{s)  and  i/o(s)  respectively  as  indicated  in  Figure 
3.1  for  3  values  of  cj.  The  inner  curve  bounds  5,,  the  subset  of  5  on  which  one  declares 
bioequivalence  and  the  outer  curve  bounds  5o,  the  remaining  part  of  S  which  leads  to 
the  declaration  of  nonequivalence.  In  between  is  the  continuation  set  C. 


c  small 


c  medium 


c  large 


Figure  1.  A  schematic  representation  of  the  optimal  boundaries  of  Problem  1,  for  various 
values  of  c. 

Since  S  is  monotone  increasing  in  c  there  are  several  critical  values,  depending  on 
which  problem  v.'e  consider,  where  S  changes  character  somewhat.  For  c  <  Cg  the  outer 
curve  does  not  close  and  (0,  s)gC  for  arbitrairy  large  values  of  s.  If  we  let  Si  and  Sg 
be  the  values  of  s  where  y,  and  yg  are  zero,  then  we  define  Ci  as  the  infinimum  of 
those  values  of  c  for  which  i,  =  Sg.  In  Problem  1  all  points  of  y^  +  s  =  1  with  the 
possible  exception  of  (0,1)  are  in  C.  In  Problem  2  the  same  holds  for  the  points  where 
=  li  with  the  possible  exception  of  (0,Tr/2). 

Hwang  tabulated  s,  and  Sg  as  functions  of  c  and  calculated  Cg  and  Ci. 
He  found,  for  Problem  1,  cn  w  0.75  and  cig  «  0.057.  Our  bounds  indicate  that  for 
0  <  Cl  <  1,  ii,  >  cj^^,  and  for  1/4  <  Ci  <  1,  sio  <  cil{2c^^  —  1).  These  imply  ci^  <  1 
and  cio  <  1/4.  Other  bounds  yield  ci,  >  (2/tc)^/^  =  0.484  and  cio  >  0.0554.  Two  of 
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these  four  bounds  are  poor,  and  one  is  good. 

Similar  results  apply  for  Problem  2  where  Hwang  calculated  C2,  ss  0.75  and  C20  ^ 
0.0518.  We  derive  S2i  >  (c2/<^o)^^^  for  0  <  C2  <  7r/4  and  S20  S  ~  1) 

2tt/21  <  C2  <  7r/4,  where  ^0  =  (27r)“*/^.  We  also  find  (tt/Sc)*/^  =  0.380  <  C2i  <  7r/4  = 
0.785  and  0.0498  <  cjo  <  27r/27  =  0.233. 

When  c  becomes  leirge  the  set  C  shrinks  but  never  completely  vanishes.  However, 
in  the  discrete  time  version  of  the  problem,  the  effective  C  shrinks  still  further  and  for 
sufficiently  large  cost  of  sampling  it  may  pay  not  to  take  observations  no  matter  what 
the  normal  prior  distribution  is.  In  particular  the  region  ne2ir  the  points  (y,s)  =  (±1,0) 
representing  very  large  values  of  t,  and  hence  of  large  observed  sample  size,  disappear 
from  C  in  the  discrete  time  version  of  the  problem. 

This  is  borne  out  both,  by  Hwzuig’s  coarse  approximations  to  the  solutions  of  the  free 
boundary  problem,  and  his  illustration  of  the  discrete  time  correction  based  on  the  solution 
of  the  continuous  time  version. 

4.  Solutions  of  the  Heat  Equation 

Ir  this  section  we  will  concentrate  mainly  on  separable  solutions  of  the  heat,  equation 
which  are  of  the  form  u  =  For  these  solutions  Hn  can  be  expressed 

in  terms  of  confluent  hypergeometric  functions  which  can  be  represented  in  various  forms. 
In  particular  we  will  present  power  series  expansions  for  odd  and  even  solutions.  Also,  for 
nonnegative  integers  n,  we  present  representations  for  H„{a)  of  the  form 

(4.1)  G„(o)  =  P„(a)<^(a)  +  Q„(q)$(q) 

where  a  =  ys~^^^,  eind  Pnio)  amd  Qn{oi)  are  polynomials  in  a  For  negative  integers 
n  we  will  express  the  H„{a)  in  terms  of  ^(q)  and 

(4.2)  J(o)  =  r  . 

Jo 

Finally  we  shall  present  a  representation  appearing  in  Goursat  (1942),  and  remark  on  a 
form  of  the  solution  useful  for  asymptotic  expmisions  for  large  s. 

In  general,  because  of  the  relation  between  the  heat  equation  and  Brownian  motion, 
one  form  of  solution  of  the  heat  equation  can  be  derived  by  eissigning  a  “heat  source”  h(y) 
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to  5  =  {(y,3)  :  s  =  0}.  Then  we  have 

(4.3)  u(y,s)  =  J  ■^4>(^—^y^^f^iw)dw  =  Eh{y  +  W^), 

where  C{W)  =  7V^(0, 1),  provided  h  has  no  serious  singularities  and  does  not  grow  too 
rapidly  as  2  — ►  oo.  We  shall  later  examine  this  form  when  h.{z)  =  2"  and  |2"|  for 
nonnegative  integer  n. 

If 

<,4.4)  u{y,s)  =  s^f‘^Hn{a)  , 

is  a  solution  of  the  heat  equation  Uyy/2  =  Ua,  then 

(4.5)  +  aH'„{a)  =  n/r„(a) 


Assuming  Hniot)  is  an  even  function  of  a  with  the  expansion 


(4.6)  Hnejoi)  — 

j=0 

with  c„o  =■  ,  we  have,  equating  coefficients  in  (4.5) 


n-2j+2  l[-f +  (J -1)] 

-  (2j)(2j  -  -  2  jli  +  (j  -  1)1 


^lnd  thus 

(4.7)  7fn,(«)  =  F(-f,5;-Y) 


where 


(4.8) 


F(^,7;x)  =  1  +  ^^  + 


7(7  +  1)21  "■  7(7  + 1)---(7  + j  -  1) 


Similarly  for  the  odd  function  of  a  with  expzinsion 


(4.9) 


H, 


>(“)  =  E 

,=o 


C  Ot 


2>+l 


with  c*o  =  1,  we  have 
(4.10) 


iH+)  +  (j-i)l  . 
2  ill +  (7 -1)1 


9 


and 


(4.11) 

In  Tables  1  and  2,  we  list  the  coefficients  for  various  integer  values  of  n  and  j.  In 
these  expansions  there  is  no  need  for  n  to  be  an  integer.  Incidentally,  it  is  easy  to  see 
that  the  derivative  of  a  solution  of  (4.5)  is  also  a  solution  for  a  lower  value  of  n.  That  is 

+  aK'(a)('1.27)  =  (n  - 

It  should  be  noted  that  if  Hn{o()  is  even  then  H'„{a)  is  odd  and  vice  versa. 

We  now  introduce  an  alternative  view  of  these  separable  solution.  Let  h{y)  =  y^/n! 
for  y>0  and  0  for  y<0  with  n  >  — 1.  Then  u(y,s)  =  s”/^G„(o)  where  G„(o) 
is  a  solution  of  (4.5)  and  is  given  by 

(4.12)  G„(a)  r  +  c)”d6  = T  (f>{a  - 

n!  J-a  Jo 

Symmetric  and  odd  solutions  of  (4.5)  are  given  by 

(4.13)  G„e(a)  =  G„(Q)  +  G„(-a) 
and 

(4.14)  G„<,(o)  =  G„(a)  -  G„(-a) 

respectively.  These  correspond  to  h{y)  =  |y|"/n!  and  /i(y)  =  sgn(y)|y|"/n!  for  all  y. 
Also  G'„(q')  =  G„_i(a)  for  n  >  0.  By  integration  we  see  that 

(4.15)  Go(a)  =  $(q) 
and  hence 

(4.16)  Goo(a)  =  2[$(q)  -  1/2] 
while 

(4.17)  Goe(a)  =  1 
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Also 


(4.18)  GiCo)  =  ^(q)  +  o$(a) 

Using  the  fact  that  G„+i(a)  =  J  Gn{o‘)da  and  Gn(o)  — >  0  as  o  — ►  — oo  it  follows  that 
for  integer  n  >  0,  G„(a)  is  of  the  form  (4.1)  where  P„(o)  is  a  polynomial  of  degree 
n  —  1  and  Qni<^)  is  a  polynomial  of  degree  n.  It  is  easy  to  see  that  Pn(a)  and  Qfi(Q) 
are  alternately  even  and  odd  polynomials  in  a  of  degree  n  —  1  and  n  respectively  with 
leading  coeflScients  1/n!.  Also  for  nonnegative  integers  m 

(4.19)  G2m,e(Q')  =  Q2m(o') 

(4.20)  G2m+l.o(Q()  =  Q2m+l(Q') 

(4.21)  =  2{P2^(o)^(a)  +  g2m(a)[$(o)  -  1/2]} 

(4.22)  G2m+l,e(Q!)  =  2{P2ni+l(“)'^(“)  +  Q2m+l(«)[^(Q()  -  1/2]} 

Thus,  for  positive  integer  n,  s"/^Qn(«)  Is  a  solution  of  the  heat  equation  which  is  a 
polynomial  in  s  and  y  of  degree  n  in  y.  For  example 

s”Qo(<i)  =  1  , 

=  y 


and 

sQ2{q)  =  s(q^  +  l)/2  =  (y^  +  s)/2 
Another  way  of  looking  at  (3n(o)  is  that 

(4.23)  Q„(q)  =  G„(q)  +  (-1)'*G„(-q)  =  ^J5(a  + 

nl 

where  £(W^)  =  iV(0, 1).  Thus,  for  integer  n>0 

[n/2] 


(4.24) 


U  ("  -  2j)'  2'j! 
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In  particular  for  n  =  2m 


(4.25) 


G'2m,e(0)  =  Q2m(0)  = 


1.3-5- ••(«-!) 


2"*m! 


Note  that  for  n  =  (2m  +  1) 


(4.26) 


(0)  =  2^(0)Pj„+,(0) 

(2m  +  1)\  Jo 


Let  V  =  12.  Then 


G2m+1,«(0)  = 


2^(0) 
(2m  +  1)! 


re-'>{2vY 

Jq 


G2m+l,e(0)  =  2^(0) 


2'"m! 


2^(0) 


(2m +  1)!  1.3-5---(2m  +  l) 


(4.27) 


i(0)  = 


^  ^  ---r.v  /  1.3-5---(2m  +  l)  ■ 

It  is  obvious  that  (7„e(a)  =  Gne(0)jFfne(o).  Similarly  G„o(a)  is  a  multiple  of  Hnoiot)- 
That  multiple  is  clearly  G'„^(0)  =  G„_i,e(0),  and  Gnoia)  =  G„_i,e(0)£r„o(a). 

The  polynomials  Pn(o()  can  be  calculated  in  a  number  of  ways.  Direct  integration  is 
one  way.  Using  the  equation  G'„(o)  =  G„-.i(a)  which  implies  that 

(4.28)  P„_i(q)^(q)  +  Q„_i(a)$(a)  =  <^(a)[P:(Q)  -  qP„(q)  +  Q„(a)]  +  Q;(q)$(q) 

gives  an  equation  for  Pn(Q;)  in  terms  of  P„_i,P^  and  Q„.  We  can  equate  coefficients 
using  the  fact  that  Pn(0)  or  P^(0)  is  known  or  that  the  leading  coefficient  is  1/n!. 
Perhaps  the  most  efficient  method  is  the  following  indirect  one  which  yields  Pn(a)  as  the 
byproduct  of  an  asymptotic  result. 

Since 

n!G„(— q)  =  /  v^<j>{v  +  o)dv 

Jo 

=  <i>{a)  J  e-’’'/2e-'“u"du 

Jo 
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we  have,  as  o  — ♦  oo 


»+l+2» 


(4.30)  G„(-»)  ~  ^ 

j=0  “ 

For  n  =  0  this  implies  the  well  known  asymptotic  nonconvergent  expansion  as  a 

(4.31)  . 

^  /V  /v3  ^  ’  92, 1 


oo 


a  a-  O'* 


2=0 


2}j\ 


Using  this  expansion,  it  follows  that  the  nonzero  coefficients  of  n!P„(a)  can  be  calculated 
as  the  first  [(n  +  l)/2]  coefficients  of  the  formal  exptinsion  of 


(4.32) 

Pniv)  ~  qniv)a{v) 

where 

(4.33) 

2=0  “ 

and 

(4.34) 

.  /„\  (2;)!  . 

The  coefiicients  of  Pn(o)  and  Qn(tt)  are  listed  in  Tables  7a  emd  7b. 

Because  Go(a)  =  $(a),  G'o{a)  =  ^(a)  is  an  even  solution  of  (4.5)  for  n  =  — 1. 
But  the  even  solution  for  n  =  0  is  the  constant  1  and  its  derivative  does  not  provide  a 
second  independent  solution  for  n  =  —  1.  However  it  is  easy  to  see  that  that  is  furnished 
by  J(a),  which  is  defined  in  (4.2),  and  is  odd.  Thus  successive  solutions,  £jtemating  even 
and  odd,  for  negative  integer  values  of  n,  are  derived  from  successive  derivatives  of  ^ 
and  J.  The  derivatives  of  4>  can  be  expressed  by 

(4.35)  ^<”*>(q)  =  Rm{oi)<t>{a) 

where  Rm{ot)  is  a  polynomial  of  degree  m  with  leading  coefficient  (  —  1)”*,  and 
corresponds  to  Hnia)  for  n  =  — (m+1).  In  fact  the  Rmiot)  are  the  Hermite  polynomials 
except  for  the  factor  (—1)”*.  Substituting  in  (4.5)  we  have 

(4.36)  R'M  -  =  -mRM 
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aaid  expanding  Rm,  for  even  m,  as 


(4.37) 
We  have 

(4.38) 


Rmi^)  —  ^  ]  QmjQ  ^ 


-[m  +  2  -  2j] 


which  implies  that  the  amj/amo  =  (—!)■' Cmj  for  even  m.  Similarly  for  odd  m, 


(4.39) 


RM  =  E 


(4.40) 


^  -[m  +  1  -  2j] 
"‘J  (2i)(2i  +  1) 


and  amj/omo  =  (-IPCj  for  odd  j  also. 

For  {Kjsitive  integers  m,  the  coefRcients  of  Rmick)  are  proportional,  with  sign  chang¬ 
ing,  to  the  coefficients  of  Qm(<^)*  Since  the  leading  coefficients  of  Rmi^t)  alternate 
between  -1-1  and  -1  and  those  of  m!Qm(o)  are  1,  it  is  clear  that  the  coefficients  of 
m!Q(a)  are  simply  the  absolute  values  of  those  of  the  Hermite  polynomials  and 


{4.41) 


m'  1 

«„(<.)  =  (-!)”■  E(-')'(;;rr2i)!  237!“’ 


We  tvirn  to  J(a),  As  a  0 


OC^  Q.^ 


(4.42) 


J(a)~a-y-^--... 


This  expansion  coincides  with  that  of  H-i^o{a)  =  aF(— 1,3/2;  — o^/2)  corresponding  to 
(4.11).  For  Q  ^  oo 

J(a)=  /  e(**-“*>/Vx=  r 

Jo  Jo 

=  rc-““+"'‘/2da 

Jo 

Jo 

-1  r'  -j  t'" 

=  n  /  e  1  H-  —z  -I-  —7  -[-... 

Jo  2q2  Sa^ 


1: 


(4,43) 


J(a)  ~  a"‘{l  +  o"*  +  3o”*  +  15a“‘  + . . .) 


From  (4.2)  it  is  easy  to  see  that 


(4.44) 


J'(ot)  =  1  —  aj(a) 


(4.45) 


J"(a)  =  -a  +  ia^ -l)J{a) 


(4.46) 


=  (o2  -  2)  -  (q^  -  3a)J(a) 


Thus 


(4.47) 


j('")(o)  =  S„(a)  +  r„(a)J(o) 


where  Sm(ot)  and  Tm(oi)  are  polynomials,  the  coefficients  of  which  evidently  relate  to 
those  of  m!Pm(o)  and  m!Qm(Q')  except  for  sign  changes.  Substituting  in  (4.5),  we  have 


(4.48) 


r;;(o)  -  ar;(a)  =  -mr„(a) 


from  which  it  follows  that  Tmioi)  =  Umioi)  because  the  leading  coefficients  coincide. 


By  formed  differentiation  of  (4.43)  we  have,  as  a  — »  oo 


(4.49) 


j('")(q)  ~  (-l)"*a-’"-^  V 

^  ^  ^  ^  ^0  (2;)!  2^j! 


Thus  the  coefficients  of  Sm(a)  are  (—I)"*  *  times  the  first  [(m  +  l)/2]  of  the  formal 
expansion  of 

tm(*^)jm(^) 


where 


(4.50) 


4^  2Jjl 
j=o  •' 


(4.51) 


[wi/2)  ,  \  (r\  "VI 
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Relating  these  expansions  to  those  of  (4.32)  -  (4.34),  we  see  why  the  coefficients  of  5^(0;) 
coincide,  except  for  sign,  with  those  of  m!PTO(o()-  The  coefficients  of  Rm  =  Tm  and  Sm 
are  given  in  Table  8. 


Comparing  Gne»  G„o,  and  with  the  confluent  hypergeometric  functions, 

we  have 

J'")(a)  =  j(”)(0)F(^,i;-y)  for  n  odd 

J<"-'>(o)  =  jl"){0)aF(2.±^,  I',—) 

d'">(a)  =  d<"-^‘'(0)aF(!^,2.,-^) 

(4.52)  d<"''(<»)  =  d*"«'(0)F(=^,i;-y) 

where 


(4.53) 


G2m,e(0) 


G2m-l,e(0) 


1 

m!2"» 

_ 2^0 _ 

1.3- 5- • -(2771  +  1) 

(-ir 

m!2”‘ 


,^(2m)(o)  ^ 


(-l)"<^o 

l-3-5----(2m-  1) 


and  (f>o  is  an  abbreviation  of  ^(0)  =  (27r) 


For  our  asymptotic  expainsions  for  large  s  in  Section  6.3  we  will  use  a  variation 
of  Equation  4.3  involving  fimctions  which  behave  like  h{y)  =  y~^  log  y^  for  large  y. 
Equation  (4.3)  corresponds  to  the  effect  at  time  s  of  a  distribution  of  heat  given  by  h. 
The  variation  we  propose, 


u(y,s)  =  s 


=  c-l/2 


^  -2<i>{-~^h{w)dw 


(4.54) 


=  J  ^[^(a^  —  a)  +  ^(a^  +  a)]  —  2^(a)^/i(a^\/s)^®^ 
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corresponds  to  a  distribution  of  heat  analagous  to  a  source  and  compensating  sink  at 
y  =  0.  This  variation  permits  us  to  deal  with  the  behavior  of  h  for  large  s  when  o' 
is  small.  Oddly  enough,  the  behavior  of  h  for  small  values  of  y  is  not  important  and  h 
can  be  modified  to  be  zero  when  |y|  <  1  with  Uttle  overall  effect.  It  is  easy  to  see  that 
for  large  s,  we  obtain 

(4.55)  u(y,5)  ~  +  . . . 

2  8  48 

and 

(4.56)  Uj,(y,s)  — /i'(s*/^Q()  +  +  ^/i^®^(s*/^(a)  +  . . .  . 

2  8 

Finally,  a  representation  of  the  solution  of  the  heat  equation  from  Goursat  (1942), 
which  was  used  in  Chemoff  (1965a)  to  justify  the  asymptotic  expansions  for  large  s,  is 
reproduced  below.  We  have 

(4.57)  u(y,s)  =  ju(yi,si)dyi  +  iu„dsi 

where  the  integral  is  taken  along  a  (yi,3i)  path  B  which  starts  at  sj  =  s  with  a  value 
of  yi  less  than  y,  passes  through  points  Si  <  s  and  terminates  at  sj  =  s  with  a  value 
of  yi  greater  than  y.  The  region  that  B  circumscribes  is  one  where  u  satisfies  the 
heat  equation.  Equation  (4.3)  may  be  regarded  as  a  limiting  case. 

4.1  Tables. 

We  present  some  tables  relevant  to  the  solutions  of  the  heat  equation.  In  Table  1  we 
present  the  even  solution 

and  in  Table  2  the  odd  solution 

(4.1.2)  = 

of 

(4.1.3)  H’:{a}  +  aH\a)  =  nHn{a) 
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as  power  series  expansions  in  a,  with  leading  coefiBcient  1.  These  are  related,  for  positive 
values  of  n,  to 

(4.1.4)  G„e(o)  =  G„(a)  +  G„{-a) 
and 

(4.1.5)  G„,(a)  =  G„(o)-G„(-a) 

These  appear  alternately  in  Tables  3  and  4.  Table  3  presents 

(4.1.6)  Q„(a)  =  G„(a)  +  (-l)"G„(-a) 
and  Table  4  gives  the  expansion  for 

(4.1.7)  V;.(a)  =  [G„(a)  +  (-l)"+*Gn(-o))/2«io. 

For  negative  values  of  n,  Hm  and  Hno  are  related  to  J^~"“^)(o)  and  ^^“"“^^(q), 
for  which  the  expansions  appear  in  Tables  5  and  6. 

The  relations 

(4.1.8)  G„(a)  =  P„(o)^(o)  +  Qniamot) 

(4.1.9)  j(")(o)  =  i2„(a)J(a)  +  S„(a) 
and 

(4.1.10)  ^(">(a)  =  P„(o)^(o)  , 

where  P„,  Q„,  P„  and  5„  are  polynomials  in  a,  lead  to  the  Tables  7  and  8  for 
Pin  Qn;  n\Pn  n!Q„,  P„  £md  S„.  Note  that  the  coefficients  of  P„  and  (?„  are  the 
absolute  values  of  those  for  P„  and  S„,  and  that  the  (— 1)"P„  are  Hennite  polynomials. 
Also  Qti(o)  is  a  solution  of  (4.1.3)  and  thus  s"^^Q„(a)  satisfies  the  heat  equation.  The 


aisymptotic  expansions,  for  a  —*  oo, 

(4.1.11) 

and 

(4.1.12) 

J(a)  ~  a“*  {1  +  o“^  +  3q“^  +  15q“®  +  . . .} 

give  rise  to  Table  9  which  presents  the  asymptotic  expansions  for  large  a  of  G„(— a)  and 
J^"^(a).  Equation  (4.1.10)  is  adequate  to  deal  with  $^"^(q)  for  large  a,  and  G„(a) 
is  described,  for  large  a,  by 

(4.1.13)  G„(a)  =  Q„(a)  -  (-l)"G„(-a) 

Finally,  in  Table  10  we  indicate  how  G„e,  G„o,  and  are  related. 
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Table  1.  The  coefficients  a^r  of  Hneio) 


n\r 

0 

1 

2 

3 

-10 

1 

-10/2 

120/24 

-1,680/720 

-9 

1 

-9/2 

99/24 

-1,287/720 

-8 

1 

-8/2 

80/24 

-960/720 

-7 

1 

-7/1 

63/24 

-693/720 

-6 

1 

-6/2 

48/24 

-480/720 

-5 

1 

-5/2 

35/24 

-315/720 

-4 

1 

-4/2 

24/24 

-192/720 

-3 

1 

-3/2 

15/24 

-105/720 

-2 

1 

-2/2 

8/24 

-48/720 

-1 

1 

-1/2 

3/24 

-15/720 

0 

1 

0 

0 

0 

1 

1 

1/2 

-1/24 

3/720 

2 

1 

2/2 

0 

0 

3 

1 

3/2 

3/24 

-3/720 

4 

1 

4/2 

8/24 

0 

5 

1 

5/2 

15/24 

15/720 

6 

1 

6/2 

24/24 

48/720 

7 

1 

7/2 

35/24 

105/720 

8 

1 

8/2 

48/24 

192/720 

9 

1 

9/2 

63/24 

315/720 

10 

1 

10/2 

80/24 

480/720 

r=sO 


4 

26,880/40,320 

19,305/40,320 

13,440/40,320 

9,009/40,320 

5,760/40,320 

3,465/40,320 

1,920/40,320 

945/40,320 

384/40,320 

105/430 

0 

-15/40,320 

0 

9/40,320 

0 

-15/40,320 

0 

105/40,320 

384/40,320 

945/40,320 

1,920/40,320 
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Table  2.  The  coefficients  a„r  of  ff„o(o)  =  . 

r=0 


n\r 

0 

1 

2 

-10 

1 

-11/6 

143/120 

-9 

1 

-10/6 

120/120 

-8 

1 

-9/6 

99/120 

-7 

1 

-8/6 

80/120 

-6 

1 

-7/6 

63/120 

-5 

1 

-6/6 

48/120 

-4 

1 

-5/6 

35/120 

-3 

1 

-4/6 

24/120 

-2 

1 

-3/6 

15/120 

-1 

1 

-2/6 

8/120 

0 

1 

-1/6 

3/120 

1 

1 

0 

0 

2 

1 

1/6 

-1/120 

3 

1 

2/6 

0 

4 

1 

3/6 

3/120 

5 

1 

4/6 

8/120 

6 

1 

5/6 

15/120 

7 

1 

6/6 

24/120 

8 

1 

7/6 

35/120 

9 

1 

8/6 

48/120 

10 

1 

9/6 

63/120 

-2,145/5,040 

36,465/362,880 

-1,680/5,040 

26,880/362,880 

-1,287/5,040 

19,305/362,880 

-960/5,040 

13,440/362,880 

-693/5,040 

9,009/362,880 

-480/5,040 

5,760/362,880 

-315/5,040 

3,465/362,880 

-192/5,040 

1,920/362,880 

-105/5,040 

945/362,880 

-48/5,040 

384/362,880 

-15/5,040 

105/362,880 

0 

0 

3/5,404 

-15/362,880 

0 

0 

-3/5,040 

9/362,880 

0 

0 

15/5,040 

-15/362,880 

48/5,040 

0 

105/5,040 

105/362,880 

192/5,040 

384/362,880 

315/5,040 

945/362,880 

20 


Table  3  G„«(a)  =  Q„(o)  =  G„(a)  +  (-l)"G„(-o)  = 


o 


o 

o 

00 

00 

CM 

to 


O^CSC0^i0«b*000>O 

e 


^  o^€^e^*^t«<or^ooo> 

s 
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420  210 

1,260 

4,725  3,150 


Tables: 

;  =  0 

Onr 


r 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

1 

-1/3 

1/15 

-1/105 

1/945 

1 

1 

-1 

1/3 

-1/15 

1/105 

2 

-2 

4/3 

-2/5 

8/105 

-2/189 

3 

-2 

4 

-2 

8/15 

-2/21 

4 

8 

-8 

16/5 

-16/21 

8/63 

5 

8 

-24 

16 

-16/3 

8/7 

6 

-48 

64 

-32 

32 

-16/9 

23 


24 


o 

o 

00 

00 


o 

00 

00 


Od 


00 


b- 


k 

« 

k 

* JB 


O  00 

^  to 

o  eo 


lO 


o 


s  s 


00 

00 

bO 


•  e 


O  c*l  Ci 

CM  «o  o> 

^  CM  M 

Zi*  ^  ^ 

CM 

eo 

CM 


8  S  8 

o  ^ 


o 

CM 


kO 

ec 


•  e 

wa 


lO 

QO 


00 


kO 


to 

CM* 


CM 


o 


O  CM 

o> 

kO 

CM 

OO  lO 

to 


00 


kO 

Oi 


CM 

kO 


o 

00 


o 

CM 

o 


CO 

CO 


CM 


kC 

Oi 

00 


o  ^ 


—  CM 
CO 

CO 

o> 


kO  kO  kO 

•M  o  ^ 

^  —  o 


o 


CM  00  00 


00 

eo 


»-<CM<0^k0t0t^000>O 


^CMeo^kOtor*ooo>o 
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mie  8.  =  Rnia)<l>ia) 

J^"\a)  =  Rn{Q)J{a)  +  Sn{ci) 


I 


•-4  ce 
CO 


^  00 
<N 


4A 


M  oO 

»  c- 

eo 


uo  o  o 

-N  40 

•  CM 


to  o 

O  CD 

—  CM 


o  tn 

CM  CM 


o>  r^ 

*  00 


^  to  CO 

•  CO 


to  to  to 


CM  00  00 

I  ^  00 

I  CO 


^  Ci’-^CMe0^i0CDr»00D>O 

e 


*-<CMC0'^t0CDt^000)O 
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279  -185  27 

-976  345 

-2,895  2,640  -588 


Table  9.  Asymptotic  Expansioos  for  large  a  for  and  C7n(— or) 


Gn{o)  ~  o' 


Onr 

n\r 

0 

1 

2 

3 

0 

1 

1 

3 

15 

1 

-1 

-3 

-15 

-105 

2 

2 

12 

90 

840 

3 

-6 

-60 

-630 

-7,560 

4 

24 

360 

5,040 

75,600 

3 

-120 

-2,520 

-45,360 

-831,600 

6 

720 

20,160 

7 

-5,040 

-181,400 

8 

40,320 

9 

-362,880 

bnr 

n\r 

0 

1 

2 

3 

0 

1 

-1 

3 

-15 

1 

1 

-3 

15 

-105 

2 

1 

-6 

45 

-420 

3 

1 

-10 

105 

-1,260 

4 

1 

-15 

210 

-3,150 

5 

1 

-21 

378 

-6,930 

6 

1 

-28 

630 

-13,860 

("+')^(o)E6„ra-’" 
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Table  10.  Relations  among  H,  J,  4t,  G  and  T 


B-2m,e  ~  *V®2(m— 1)  ^2m,e  ~  |®2in(^Jm,e  —  |®2m|Q2m 

H-im-Xfi  =  lO’im  Bim+lfi  =  W2m\G2mHfi  =  |02m|Q2m+l 

B-2m-l,e  =  ^2m-l,e  =  jfr2m|C^2m-l,e20O  =  |i2m|72m-l 

R-imfi  =  ^^^”*~*V^2m^  Bxmfi  —  l^m|G’2m-l,e2^  =  |b2m|72m 


n  0  2  4  6  8  10 

a„  1  -2  8  -48  384  -3,840 

bn  1  -1  3  -15  105  -945 


aim  =  (-ir2"‘m!  62m  =  (-1)”*^  =  •  3  •  •  •  (2m  -  1)) 
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5.  Bounds 


In  this  section  we  derive  some  qualitative  properties  of  the  optimal  sequential  rules 
for  Problems  1  and  2.  Most  of  these  establish  the  bounds  described  in  Section  3.  First  we 
have 

Proposition  1.  The  optimal  continuation  set  C  decreases  with  cj  in  Problem 
1,  and  with  in  Problem  2. 

Proof:  Let  di2,/)  and  C  correspond  to  ci  and  and  C*  to  cj  >  c  and  let 

(y,s)€C*.  Then  there  is  a  stopping  rule  S<s  for  which 

dl2{y,s)>E{d\2[yiS),S]\Y{s)  =  y} 

=  E{du[YiS),S]  +  (ct  -  c,)S-Mr(s)  =  y) 

>  piy,s)  +  {c*  -ci)s~^ 

and 


di2{y,s)  =  dl^iy.s)  -  (cj  -  Ci)s"^  >  p{y,s) 


Hence,  (y,s)  €  C  and  C  DC*. 


The  scime  proof  applies  for  Problem  2. 


We  state  two  related  propositions  and  a  more  general  lemma,  on  which  the  proofs  of 
the  propositions  are  based. 

Proposition  2.  In  Problem  1  with  stopping  cost  di2,  the  optimal  continuation 
set  Cl  contains  {(y,s) :  y^  +  s  =  1,  s  >  0,  y  ^  0}. 

Proposition  2*.  In  Problem  2  with  stopping  cost  d22,  the  optimal  continuation 
set  C2  contains  {(y,s) :  =  1,  s  >  0,  y  ^  0}. 

Suppose  that  yo(s)  is  a  differentiable  curve  along  which  d{y,s)  is  continuous  and 
differentiable,  but  along  which  the  right  hand  and  left  hand  derivatives  with  respect  to 
y,  dy  and  d~  are  defined  finite  but  unequal. 
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Lemma  1  For  C{W)  =  N{Q,  1),  E{d[yo{s)  +  WVI,  s  -  6]  =  d(yo{s),  5)  +  - 

d;)VS  +  oiy/S). 


Proof; 

dlyo(s)  +  Wy/6,s  -  6]  =  d(yo(s  —  6)  +  W\/6  +  yo{a)  -  yo(s  -6),s  ~  6] 

=  d[yo(5),  s]  +  [wV6]d^  +  OiS)  for  W  >  0(y/S) 

=  d[yo(s),  5]  +  [W\/S]d;  +  0{6)  for  W  <  0(y/S) 

Hence 

F{d[yo(s)  +  Wy/Sy  a  -  ^]}  =  d[yo(s),  a]  +  y/S{d^  f  w4>{w)dw  +  d~  f  w<l>{w)dw}  +  0{S) 

Jo  J-00 

=  d[yo(s),  5]  +  <f>oy/6[di  -  d;]  +  0{6) 


Note:  Implicit  in  the  proof  above  was  the  assumption,  not  clearly  stated  here,  that  d  does 
not  grow  too  fast  as  y  goes  to  ±00. 

Propositions  2  and  2*  are  corollaries  of  Lemma  1.  For  Proposition  2,  ^2y  =  0,  and 
d'[2y  =  2y  for  y^  +  s  =  1  and  y  >  0.  While  df2j,  =  0  and  =  2y  for  y^  +  s  =  1 
and  y  <  0  .  In  both  cases  dj’jj  —  dfjy  <  0  and  the  procedure  of  letting  S  =  s  — 6  leads 
to  a  reduction  of  risk  over  stopping  at  (y,5). 

For  Proposition  2*  we  have  for  s*/^Gie(ys”^^^)  =  1  and  y  >  0,  ^ 

d^2y  =  2$(ys”^/^)  —  1  >  0  For  y  <  0,  d22j,  =  0  and  d^2y  —  2$(ys“*/*)  —  1  <  0.  Once 
more  we  have  <^2y  ~  <^22y  <  0  in  each  case  and  the  result  follows. 

■ 

Propositions  2  and  2*  fail  to  determine  whether  the  points  (0, Si)  =  (0,1)  and 
(0,32)  =  (0, 7r/2),  where  sy^G\^{Q)  =  1,  are,  or  are  not,  in  the  continuation  sets.  This 
question  is  partially  addressed  by  Propositions  3  and  3*. 

Proposition  3.  If  ci  <  (2/n-e)^/^  =  0.484,  (0,1)  e  C  in  Problem  1. 

Proposition  3*  If  C2  <  (tt/Sc)^/^  =  0.380,  (0, 7r/2)  €  C  in  Problem  2. 
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Proof:  With  £(W)  =  iV(0, 1),  we  have 


/OO 

min(l,u)^6  +  1  —  6)4{w)dw 

-OO 

=  Cl  +  1  +  ^{ci  +  J  (w^  -  l)<^(w)clw}  +  o(S) 


=  d(0, 1)  +  ^{ci  —  2^(1)}  +  o(S) 

Hence  (0,1)€C  if  ci  <2^(1). 

d22(fyy/s,  I  -  <)  = 

Since  Gie(Q)  =  2^o(l  +  +  o(o^)  for  ot  — »  0 

E^d22{WVS,  f  ~  ^)}  =  ^*=2  +  +  o(^) 

=  <f(0,7r/2)  +  ^1^02  +  J  ^  +o(^) 

Thus  (0, 7r/2)6C  if 

2  1  fl 

C2  <  -^ —  /  (w^  —  l)^{w)dw  =  ^^(1)  =  -s/ir/^t 
4  TT  y_i  2 


From  a  slight  perturbation  of  Propositions  3  and  3*  it  follows  that  ci,-  >  (2/7rc)*/^  and 
C2j  >  (tt/Sc)^/^.  Upper  bounds  on  ci^  2uid  C2i  follow  from  Propositions  4  and  4*. 


Proposition  4.  In  Problem  1  all  points  (0,s)  with  0  <  s  <  are  in  S  if 
0  <  Cl  <  1.  If  Cl  >  1,  (0,s)  €  5  for  all  s  >  0,  and  ci,-  <  1. 


Proposition  4*.  In  Problem  2  all  points  (0,  s)  with  0  <  s  <  (02/^0)^^^  are  in 
S  if  0  <  C2  <  ?r/4.  If  C2  >  7r/4,  (0,s)g5  for  all  s  >  0,  and  hence  C2i  <  ?r/4- 

Proof:  For  Problem  1,  let  0  <  si  <  s  <  cj^^  <  1.  Then 

dl2(0,s)  =  Ci5~^  +  s 


but  for  all  y 


<ii2(y,si)  >  CiSf^  +si  >  di2(0,s) 
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since  Cis  ^  +  s  is  decreasing  in  s  for  0  <  s  <  c|^^.  Hence  it  is  diszulvantageoiis  to 
continue  from  (0,s). 

K  Cl  >  1  and  0  <  si  <  s,  then  for  all  y 

dniy^st)  >  di2(0,si)  =  cisf*  +inin(l,si) 

which  is  a  decreasing  function  of  Si.  Hence  di2(y,Si)  >  <ii2(0, s)  and  (0,5)  e  5.  It 
follows  immediately  that  ci,  <  1. 

For  Problem  2,  let  0  <  52  <  5  <  (27rc2)*/^  <  7r/2.  Then 

<f22(0.^)  =  C25“*  +  v^2/7r5*/^ 

since  Gie(a()  assumes  its  minimum  value  of  2^o  =  \/2/fl’  at  ot  =  0  and  2<f>os^^^  <  1. 
For  all  y 

d22iy,S2)  >  C25^*  +  >  d22(0,s) 

since  C2s“^  +  is  decreasing  in  5  for  s  <  (27rc|)^/®.  Once  again  it  is 

disadvantageous  to  continue  from  (0,5). 

If  C2  >  ir/4  and  0  <  S2  <  5,  then  for  all  y 

d22iy,S2)  >  ^22(0,52)  =  C257*  +  min(l,  y/^sy^) 

which  is  a  decreasing  fimction  of  52.  Thus  ^22(1/,  52)  >  <^22(0, 5)  and  (0,s)  G  S,  and 
consequently  62, •  <  7r/4. 

■ 

Some  explicit  bounds  for  Sio  and  520  are  given  in  Propositions  5  and  5*. 

Proposition  5.  For  1/4  <  ci  <  1,  (y,5)  G  5  if  5  >  ci/(2c|^^  —  1). 

Proposition  5*.  For  2-k/21  <  C2  <  rr/4,  (y,s)  G  5  if 

5  >  C2/(3cJ/VJ^"  -  1) 
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Proof:  Let  5  be  an  arbitrary  stopping  procedure  (starting  from  (y,s))  for  Problon  1. 
Let  A)  represent  the  contribution  of  the  event  A  to  EX,  \.c.,  Jy^X(w)dP{w). 

Then 

=  i;{dx2[n5),  S]}  =  E{di2[Y{S),  S);  S  >  1}  +  E{dn[Y{S),  S];  S  <  1} 

>  (cis-^  +  1)P{S  >  1}  +  2cJ/^P{S  <  1} 

since  for  s  >  S  >  1,  di2(y,S)  =  ciS“*  +1,  and  for  0  <  5  <  1,  di2(y,S)  >  cjS”*  +  S  > 
2c;/^ 

As  Cl  ranges  from  1/4  to  1,  ci/{2cj^^  -  1)  decreases  monotonically  from  oo  to 
1  and  hence  s  >  1,  and  di2(y,s)  =  cis~^  +  1,  and 

Pi  -  d^2iy,  s)  >  i2cY^  -  c,s-»  -  1)P{S  <  1}  >  0, 

and  hence  (y,  s)  G  S. 

Let  S  be  an  arbitrary  stopping  procedure  for  Problem  2.  Then 

P2  =  P{d22[y(5),5]}  =  P{d22[r(5),5);S  >  7r/2}  +P{f/22[>'(P),*?];P  <  ’r/2} 

>  (C2S-^  +  1)P{S  >  ir/2}  +  (27c2/2ir)*/^P{S  <  7r/2} 

since  for  s  >  S  >  vj2,  d22(y,5)  >  C25"^  +  1  and  for  0  <  S  <  7r/2,  d22(y,  5)  > 
C2S-1  +  5^/2Gie(0)  =  C25-1  +  2«ioS‘/2  >  (27c2/27r)>/’. 

As  C2  ranges  from  27r/27  to  7r/4,  decreases  monotonically  from 

00  to  7r/2  and  hence  s  >  7r/2  and  d22(yi'S)  =  C2S~’ +  1,  and 

P2  -  d22{y,s)  >  [(27c2/27r)»/^  -  cj^"^  -  1]P{S  <  7r/2}  >  0 

and  hence  (y,  s)  6  S. 


1  /2 

As  a  consequence  of  Propositions  5  and  5*,  it  follows  that  iio  ^  Ci/(2cj'  —  1)  and 
520  <  C2/(3c5/V?^’  -  1),  CM  <  1,  C2i  <  7r/4,  c,„  <  1/4,  and  C2„  <  27r/27. 

In  Propositions  6  and  6*  we  derive  inner  bounds  on  C  for  Problems  1  and  2  for  ci 
and  C2  suflSciently  smzdl.  These  incidentally  imply  that 

5.„>c.,=  {4_^m|cja^  +  l)[l-^)}  =0.0554 
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where 


J'(qo)=  inf  J'(a)  = -0.2847 

0<Of<oo 

and 


ao  =  2.124. 


Also 


=  C?,(a)|l  -  J'(o)/.''(<>o)))  =0.0498 

2/  ^0<a<ao  J 


Proposition  6.  For  0  <  ci  <  C13  in  Problem  1,  the  optimal  continuation  set 
contains  Ai  which  consists  of  the  set  of  all  (y,  s)  for  which  |a|  <  oq  and 

1  +  {1  -  4c„(l  -  -  r(a)/ J-(ao)l)‘/^ 

2(1  +  a2) 


s  > 


Proposition  6*.  For  0  <  cj  <  £33  in  Problem  2,  the  optimal  continuation  set 
contains  As  which  consists  of  all  (y,5)  for  which  |q|  <  ao  and  s  >  ^^.(a) 
where  a+(a)  is  the  larger  of  the  two  positive  roots  of 

,  +  C23(l  -  J\a)/J'iao)]  =  0  . 


Proof.  Let  u(y,  s)  =  1  +  cis“V'(o)/J'(qo),  the  second  term  of  which  imitates  Cis“* 
when  s  -+  0  and  a  =  qq.  Our  proof  involves  showing  that  u  <  di2  on  a  set  A,  and 
u  =  di2  on  the  boimdary  of  A.  Thus  u  is  the  Bayes  Risk  corresponding  to  the  procedure 
defined  by  C  =  A.  Since  p{y,s)  <  u(y,s)  everywhere,  then  p  <  d\2  whenever  u  <  di2 
and  hence  the  optimal  continuation  set  Cj  D  A.  Now  Ci  is  monotone  decreasing  in  cj . 
But  our  set  A  is  monotone  increasing  in  cj.  Hence  Ci  must  contain  Ai  =  A(ci3). 

Since  all  points  (0,  s)  for  which  s>l  are  included  in  Cj,  it  follows  that  cio  >  C13. 
Figure  2  helps  to  illustrate  the  details  of  the  proof,  which  follows. 

The  two  curves  y^  =  OqS  and  y^  +  s  =  1  intersect  at  (±yo,5o)  where  so  = 
{oq  +  1)"^  <  1  and  7/0  =  The  set  Bi  =  {(y,s) :  y^  +  s  >  1,  y^  <  Oqs)  will  be 

part  of  A.  The  two  points  (±yc,so)  will  be  connected  by  a  curve  segment  s+(y)  which 
together  with  s-(y)  sxe  the  roots  of 

Q{s)  =  (a^  +  l)s^  -  s  +  cx[l  -  J'ia)/J'{ao)]  =  0. 
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These  roots  are  real  when  the  discriminant 


D  =  l-  '''^^  +  l)cill-ria)/J\ao)]>0 

for  |a|  <  Qo.  That  is  the  case  because  0  <  ci  <  c\3.  Moreover  0  <  D  <  \  and  the  roots 
5+  and  3_  are  given  by  (1  ±Z>^^^)/2(1  +a^)  satisfying  0  <  s_  <  <  (o^  + 1)“^  <  1. 

Let  A  =  BiUB2  where  =  {(y,s)  :  <  OoS,s  >  s'*"}. 


Figure  2.  Regions  bounded  by  |qo|  =  a,  «  =  di2,  and  Q  =  0  for  Problem  1  for  3 
values  of  Ci. 

Now  we  show  that  <ii2  —  u  >  0  on  A  with  equality  on  the  boundeu-y.  On  Bi 

di2  —  u  =  cis~*[l  —  J'(a)/J'(ao)]  >  0 


and  on  B2,  where  s  > 

di2  —  u  =  s”*Qfs)  >  0. 

Equality  on  the  boundary  of  A  is  apparent.  Thus,  p(y,  s)  <  u(r,  s)  <  d(y,  s)  on  A,  Ci  D 
A.  As  Cl  increases  subject  to  0  <  ci  <  ci3,  D  and  s+  de-rease,  enlarging  B2  and 
A  which  approaches  Ai  =  A(ci3)  but  C\  decreases.  Hence  Ci  D  Ai  =  A(ci3)  D  A(ci). 
Proposition  6  follows. 


Oiir  proof  has  incidentally  shown  that  Ci  includes  all  points  (0,  s)  for  which  s  >  1 
if  Cl  <  Ci3  and  hence  Cio  >  cia. 

The  proof  of  Proposition  6*  proceeds  along  similar  lines.  We  use  the  same  solution 
u  of  the  heat  equation  with  ci  replaced  by  cj.  The  two  curves  =  UqS  and 
s^/^Gie(ys~^/^)  =  1  intersect  at  (±yo,s5)  where  0  <  sj  =  [Gie(Qj)]“^  <  7r/2  and 
Vo  =  ooSq  '  .  For  the  time  being,  let  us  assume  that  these  two  points  are  connected  by  a 
curve  on  which 


(?*(s,  a)  =  3^/^Gi,(a)  -  s  +  C2[l  -  J'(a)/7'(ao)]  =  0 


where  <  1,  and  Q*(s)  >0  if  s  >  s^/^Gie(Q)  <  1  and  |a|  <  qq. 

Then  we  let  .4*  =  Bj  U  where  Bj  =  {(y,s)  :  s*/^Gie(ys~^/^)  >  1,  <  Oqs} 

and  Bj  =  {(y,^)  :  s^^^Guiys'^^^)  <  1,  y^  <  a^s,  s  >  s^}.  On  BJ,  6.22  -  u  = 
C2S“^[1  —  J'(a)/J'(ao)]  >  0  and  on  B|,d22  —  u  =  s"^Q*(s,o)  >  0.  The  rest  of  the  proof 
follows  as  in  that  of  Proposition  6,  provided  s!j.  is  nonincreasing  in  C2.  It  remains  only 
to  demonstrate  the  existence  of  the  desired  s^. 

For  each  fixed  a  in  (-ao,oo),  Q*(si,a)  >  0  at  sj  =  1/Gie(a).  Also 

for  s  >  s5  =  4/9Gfe(a).  But 


C?*K,a)  = 


-4 

21GUa) 


+  C2 


J'jcx)  ' 
J'icio). 


<  0 


for  C2  <  C23  by  definition  of  C23.  Thus  there  is  a  value  s!j.(Q:)  =  sup{s  :  Q*(s,a)  = 
0,  5$  <  s  <  sj}.  Since  dQ^/ds  >  0  and  dQ*/da  =  _  1/2]  _  C2  J"(a)/J'(ao), 

for  s  >  Sj  and  |a|  <  Oq,  defines  a  curve  and  one  for  which  slj.  is  decreasing  in  C2. 


6.  Asymptotic  Expansions 

This  section  has  4  subsections  devoted  to  eisymptotic  expansions  of  the  optimal  boimd- 
aries  and  risks  for  Problems  1  and  2  when  s  is  small  and  when  s  is  large.  Most  of  our 
derivations  will  be  formal,  since  the  argument  required  for  rigor  are  complex  and  have  been 
outlined  elsewhere  for  similar  problems,  e.g.,  in  Chernoff  (1965). 
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6.1  Problem  1,  s  small. 

The  boundzuy  of  the  optimal  stopping  set  is  near  y  =  ±1.  To  simplify  the  calculations 
we  will  consider  solutions  of  the  heat  equation  of  the  form  u(y  —  1,5)  +  u(— y  —  1,  s)  where 
u(— y  — 1,5)  will  have  a  negligible  effect  near  y  =  1.  In  effect,  for  s  close  to  0  we  have 
transformed  the  problem  to  that  of  dealing  with  ^13(2,5)  near  2  =  0  where  2  =  y  — 1 
and  di3(z,s)  =  di2iy,s),  i.e., 

^13(2, 5)  =  Ci5"^  +  min(l,  2^  +  22  +  5  +  1) 

(6.1.1)  =  Ci5“^  +  1  +  min[0,s(^^  +  1)  + 

and  =  25“^^^.  For  further  convenience  we  subtract  several  solutions  of  the  heat  equation 
from  di3.  First  we  subtract  cis~^J'{/3)  to  reduce  considerably  the  unpleasant  effect  of 
the  5”^  singularity  when  s  and  0  are  small.  To  make  the  remaining  term  more 
symmetric  in  2,  we  subtract  1  +  +  s(/9^  +  l)/2  and  we  have 

(6.1.2)  ^13(2, 5)  =  Cis~^  d'{0)  +  1  +  ^  +  l)/2  +  cidn{z,  s) 

and 

(6.1.3)  du(^,  s)  =  5-'  [1  -  /'(/?)]  -  cf'  +  5(^2  +  i)/2| 


We  will  seek  approximations  to  inner  and  outer  boundaries  ;S,  and  0oi  that  we 
assume  are  close  to  y^  +  5  =  1  in  the  original  vsiriables,  and  a  solution  of  the  heat 
equation  u  so  that  u  and  du/d0  match  di^  and  di^p  =  ddn/d^  on  both  boimdaries. 
With  an  expansion  for  the  boundaries  of  the  form 


(6.1.4) 


/?  = 


=1/2 


■{1  +  OjS  +  025^  +  035^  +  .  .  .} 


we  can  expand  dn  and  dwp  and  appropriate  versions  of  u  and  up  in  powers  of  s. 
Letting 

V  =  sgn[s^/^j6l  +  s(l  +  /9^)/2] 


we  have 
(6.1.5) 


3  15  105  ■ 


5'/2i9  +  5(^2^1)/2 
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(6.1.6) 


=  s  ^ 


—  Cj  +  s0 


Note  that  for  small 

(6.1.7)  +  +  l)/2  =  s^ 

and  hence, 

(6.1.8)  V  =  sgn(l  —  4ai ) 
Now  we  have 


1  -4oi' 

-1. 

Oi  —  2o2 

1 

Oj  +  2o2  —  403' 

8 

T  S 

4 

"T  5 

8 

+ 


1 

Oi  1 

2 

=  i+^ 

.  2  48 

+  s^ 

(6.1.9) 


+  s^ 


Cj  Cj  Oi  1  _i  /oi  iM 

03  QiQ2  02  ^  2  ,  Ql  1  _1  /  202  —  Ol  Nl 

2  2  12  8^*  "^  160  26880  "^^^  4  )\ 


+ 


(6.1.10) 


du^  =  -s"^/^|l  +  s 

+  s’  [».  -  ja,  -  \a\  +  i„,  -  ^  -  cr'^lj  } 


-  6+cr^v 


+  S-' 


1  1 
■  2°'  80 


At  this  point  we  introduce  a  first  approximation  to  the  optimal  solution  of  the  heat 
equation  which  will  match  dn  and  on  the  inner  and  outer  boundaries  determined 

by  appropriate  terms  (oi,,02,, . . .)  and  (oi„,02o,...)  with  oi,  >  1/4  >  Oio  so  that 
Vi  =  sgn(l  -  4oii)  <  0  and  =  sgn(l  -  4aio)  >  0.  Let 

(6.1.11)  Ui  =  ris~*/^7(/3)  +  r2 

with 


(6.1.12) 


«i3  =  ris 


Expanding,  we  have 


(6.1.13) 


1  ■ 

■-01  1  ■ 

Ui  = 

r2  -  2r, 

+  ris 

.2  "^24, 

+  ris^ 

+  ri5^ 


-03  ,  ^  ,  £i  _  ^  ,  1 

2  8  8  96  13440 


1 

480 
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(6.1.14)  =  r..-A|l  -  ,  .  i  +  [-|i  +  1]  +  ,3  ^  i  -  jL]  + 

Matching  coefficients  alternately  with  and  we  have 


(6.1.15) 


ri  =  -1 


(6.1.16) 


rj  =  -1/4 


We  may  also  proceed  further  after  noting  that 


(6.1.17) 


[^j + +“■  (s + v)  “  (9S + V)] 

I  +  3’ +  <.2  +  ^)  -  a.  (2^  +  ^)  +  26^]  + 


(6.1.18) 


du0-ui0^-s  +  i  +  +  ^2  ^  _  £1^ j 

+  »’[“3-T+-(-^- V)  +  2^]+- ■•} 


By  adjoining  to  ui  the  solution  (of  the  heat  equation) 


(6.1.19) 


«2  =  +  r4s(l  +  0^) 


(6.1.20) 


=  ^  2r4S0 


we  match  the  coefficients  of  and  s  in  di40  —  ui^  —  uj/j  and  du  —  u\  —  U2 


We  have 


(6.1.21) 


(6.1.22) 


^3  ,  ,2  °1^3  , 

U2  =  s  Y  +  r4  +s  -Y-  + 


^4  ,  3  ‘*2*'3 

7J+’  l~ 


+  -^  +••• 


«2/J  =  -■» 


|sr3  +s^r4 


+  s  air4  + 
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This  leads  to 


(6.1.23)  ’■3=01  +  ^  +  Cj~‘t; 

and 


(6-1-24)  T+^‘  =  ^ 

But  the  first  equation  represents  two  equations  depending  on  v,  the  sign  of  (1  —  4ai). 
It  stwds  for 


(6.1.25) 


=  ““+12 


(6.1.26)  ^3  =  oio  +  ^  +  c~^ 

or 


(6.1.27) 


aio  =  au  -  2ci  ^ 


At  this  point  we  have  expressed  Uio,  r3  and 


(6.1.28) 


J_  _  li  =  ~°i«  _  JL  +  £l 

48  2  2  48  +2 


in  terms  of  an  undetermined  parameter  oi,  >  1/4.  Now  we  introduce 


(6.1.29)  «3  =  ^1  +  2^^  + 

with 


(6.1.30) 


U30  =  -r5S^^^(l  +  0^)  +  4r6S^ 


Expanding  along  the  boundaries  we  have 


(6.1.31) 

and 


2 

1  3 

airs 

»*6 

U3  =  5** 

[t+^ 

+  S^ 

2 

■*■24  + 

2. 

+  . . . 


(6.1.32) 


U30 


=  -.,-1/2 


{ 


rss^  +  s® 


Y+2r6 

4 
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Setting  the  coeflBcients  of  5^  amd  in  du  —  ui  —  U2  —  U3  and  di4^  — «i^  — U2j8~“3^ 
equal  to  zero,  we  have 


(6.1.33) 


V960  8  J  V  2 


2 


+  r6 


and 

(6.1.34) 


Considering  that  ai  and  v  can  take  on  two  distinct  values  in  (6.1.33)  and  r3,r4,r5  and 
rg  are  restricted  to  single  values,  it  follows  that 


(6.1.35) 


Since  Oii  —  Oio  =  2ci  ^  and  (aj,  —  Cj^)  =  (a^  +  aio)2ci  ^  =  4(ai,  —  cf’)cj  \  it  follows 
that 


(6.1.36) 
and 

(6.1.37) 


ra  =  1/3  , 


r4  =  - 


T_ 

48  ’ 


(6.1.38)  oii  =  ^  +  Cl  \ 
and 

(6.1.39)  a,„  =  i-cr'  . 
FVom  (6.1.34)  it  follows  that 

(6.1.40)  a2i  —  020  =  — cf  *  . 


Now  we  introduce 

(6.1.41)  U4  =  +  ^0^)  +  r8S^(l  +  +  ^0^) 

O  lo  10 
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and 

(6.1.42)  U4/J  =  -r7S®/^(l  +  2/9^  +  \p*)  +  r85^{6^  +  4;?’  +  . 

6  0 

Expanding  along  the  boundaries  we  have 

(6.1.43)  U4  =  s^(j+r8)  +  ... 
and 


(6.1.44) 


U4/J  =  -s  +  ...}  . 


Matching  coefficients  of  and  in  di4  —  ui  —  uj  —  U3  —  U4  and  its  derivative  with 
respect  to  we  have 


Since  the  rj  do  not  depend  on  whether  we  use  Oji  or  ajg,  the  equation  (6.1.45)  implies, 
with  some  algebraic  labor, 


(6.1.47) 


rs  = 


80 


which  combined  with  (6.1.33)  and  (6.1.34)  yield 


(6.1.48) 


2 


(6.1.49) 
and 

(6.1.50) 


167  c~^ 

960  4  ■ 
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Applying  (6.1.46)  to  aji  and  Cjo,  it  follows  that 


(6.1.51) 


and  ry  and  rg  can  be  evaluated  in  terms  of  agg.  Thus  far  we  have  for  the  inner  and 
outer  boimdaries 

(6.1.52)  i,  =  y.  -  1  =  -1  jn-  Q  +  . .  .| 

(6.1.53)  io  =  yo  -  1  =  ^  (i  ~  (4  +  •  -  I 

which  compare  with  y^  +  s  =  1  where 

(6.1.54)  y  =  (1  -  ^)-'/'  =  ^ 


The  approximation  to  the  optimal  risk  given  by  u  =  uj  +  U2  +  “3  +  . . .  is 

u  =  -l-lz-  ^{s  +  z^) 

(6.1,55)  _  +  2..^  +  1.)  +  , . . 


Deriving  further  terms  becomes  straightforward,  but  tedious,  and  requires  considerable 
care. 


6.2  Problem  2,  s  small 

In  Problem  2,  the  boundaury  is  ailso  neair  y  =  ±1  when  s  is  smaill,  and  once  again 
we  translate  to  z  =  y  —  1.  We  have  for  laurge  |q| 

GuM  =  2[^(a)  -  |a|(l  -  ^(|a|)))  +  |a| 

(6.2.1)  «|«|  +  ,«a)|^-^  +  ...} 

Thus  for  smaill  s  auid  y  close  to  one,  Gie(a)  is  very  close  to  |al  and 

d22iy,s)  =  +rain(l,s*/^Gie(o)] 

(6.2.2)  «  C2s“*  +  min(l,y)  =  1  +  i(/2  +  d23(-2,s) 
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where 


(6.2.3)  d2z{z,s)  =  CiS-^-\z\/2. 

With  the  transformation 

(6.2.4)  2*  =  02,  s*  =  c?s,  a  =  (2c2)“^^® 
we  have 

(6.2.5)  J2,(z,s)  =  (c,/4)‘/’l<i,4(2-,3')  + 
where 

(6.2.6)  d24(z,s)  =  s-‘[l  -  J'(/3)]  - 

is  sjrmmetric  in  j3,  and  therefore  easier  to  handle  than  <ii4. 

We  attack  the  asymptotic  expansion  for  ^24  near  5  =  0  and  2  =  0  by  assuming 

^  +  . . . 


in  which  case  we  can  expand,  for  positive  0 

^24  =  5^(ai  -  oi)  +  5 


<*1 


20102  -  y  -  02 


(6.2.8) 


+  s 


2O1O3  +  O2  —  2O1O2  +  ^5°!  ~  °3  I  + 


^2A0 


=  «l/2 


{ 


(2oi  —  1)  +  s' 


(2a,-ia5) 


(6.2.9) 

Let 


+  •S®  ^2o3  —  40]02  +  +  . .  .|^ 


ui  =  ris^fr4e(/3)  =  ri  5^(1  +  2/3^  +  /3^/3} 

=  ris^|l  +  s^(2o^)  +  s®(4oi02  +  Oi/3)  +  s®(4oi03 

(6.2.10)  +  2o|  4*  40J02/3)  +  . .  .^ 
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=  ris^ 

(6.2.11)  =  ris^/^|4ai  +  5^[4a2  +  4af/3]  +  5®[4a3  +40^02]  +  • 

To  match  coefficients  we  need  2ai  —1  =  0  or 


} 


(6.2.12) 

and  Oj  —  Oj  =  rj  or 
(6.2.13) 

But  then 


ai  =  1/2 


r,  =  -1/4. 


(6.2.14) 


<^24  —  Ui  =  S®  —  +5 


2  1  1  02  1 

^  6  ^  960  2  192 


+ 


+  3' 


(6.2.15)  d240  —  tti0  =  ^202  ~  g  +  9 

which  suggest  202  —  1/6  +  1/2  =  0  or 

(6.2.16)  02  =  -1/6  . 

Since  the  coefficient  of  s®  in  ^24  ~  “1  is  5/48  we  take 


1  «i 

+  80  +  T 


+ 


(6.2.17) 

2Lnd 


u.  =  ^ 

1 5a? 


l+5«^  +  |^‘  +  .. 


=  - 1-  s  — ^  -h  ...  —  s - f-  S - h 

48  48  48  192 


(6.2.18) 


=  s 


13  25 

Now  the  coefficient  of  in  ^24^  —  is  203  +  =  —  and 


(6.2.19) 


03  =  7/30  . 
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We  can  continue  in  this  manner,  taking  advantage  of  the  symmetry  of  the  problem.  In  the 
meantime  we  have  obtained 


(6.2.20) 


for  the  inner  and  outer  boimdaries  near  (y,s)  =  (1,0).  The  risk  fimction  can  be  approx¬ 
imated  by  ui-fU2.  Note  that  this  expansion  applies  to  <£34  where  the  variables  2  and 
s  refer  to  az  and  a^s  in  the  original  coordinates  with  a  =  (2c2)“'^®,  i.e., 


6.3  Problem  1,  large  s. 

For  c  large  enough,  the  continuation  region  is  bounded.  But  for  c  <  cio  it  becomes 
imbounded  and  it  is  desirable  to  find  an  asymptotic  expansion  for  the  boundary  and  risk 
ftmction  when  s  is  large  (which  corresponds  to  a  great  deal  of  prior  uncertainty  about 
/i).  Prom  Proposition  6  in  Section  5  we  know  that  for  c  <  cjs,  and  s  >  1,  that  the 
optimal  boimdary  will  have  |d|  >  oq. 

This  problem,  where 


d\2{yis)  =  cis  '-1-1  for  5>1 

represents  a  somewhat  unusual  Ccise  for  large  s.  For  a  stopping  risk  cjs”*  -hi  for  all 
s,  the  optimal  strategy  is  obviously  that  of  never  continuing.  It  is  only  the  behavior  for 
s  <  1  that  leads  to  the  possibility  of  continuing  when  s  is  large. 

In  contrast  with  this  problem,  consider  the  canonical  version  of  that  related  to  deciding 
the  sign  of  the  mean  drift  of  a  Wiener  process,  the  discrete  time  version  of  which  is  to 
decide  the  sign  of  the  mean  of  a  normal  distribution.  In  that  problem  where 

the  main  component  of  the  optimal  Bayes  risk  for  large  s  and  moderate  a  is  a  multiple 
of  There  is  another  component  which  is  relatively  small  for  Isu-ge  s  and 

moderate  a  =  But  for  large  a,  where  the  boundary  is  located,  the  first  term 

becomes  negligible  and  it  is  the  second  term  which  controls  the  location  of  the  boimdary. 
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This  fact  has  statistical  significance.  On  one  hand  the  solution  of  the  sequential 
problem  is  related  to  the  backward  induction  of  dynamic  programming,  and  depends  on 
the  structure  of  the  losses  when  s  is  small  or,  after  much  samphng,  has  reduced  the 
imcertainty  about  the  unknown  mean  fi.  The  nature  of  the  losses  after  considerble 
sampling  has  an  effect  on  our  expected  loss  and  the  optimal  stopping  rule.  On  the  other 
hand  it  has  only  a  minor  effect  on  an  early  decision  of  whether  evidence  is  overwhelming 
enough  to  stop.  In  brief,  in  statistical  decision  problem,  early  decisions  are  typically 
insensitive  to  the  consequences  of  potential  late  decisions. 

Our  current  problem  lacks  this  property.  One  consequence  is  that  our  asymptotic 
expansions,  after  the  mzun  terms,  will  depend  on  some  constant  which  depends  on  the 
nature  of  the  stopping  risk  for  0  <  s  <  1.  That  constant  is  not  available  now.  In  the 
problem  of  deciding  the  sign  of  the  mean,  bounds  on  the  corresponding  constant  were 
derived  by  Bickel  and  Yahav  (1967)  and  Mallik  and  Yao  (1984). 

First  we  drop  the  1  from  di2  and  so  we  have  dis  =  cjs"*.  Now  let  us  represent  a 
solution  of  the  heat  equation  with  the  form 

/oo 

4>{w)h{y  +  ws^^^)dw  =  roVo  +  ^ 

•OO 

while  equations  (4.54-4.56)  would  be  more  rigorous.  Here,  our  first  approximation  to  h{y), 
giving  rise  to  vj ,  is 

(6.3.2)  /ii(y)  =  riy"^logy^ 


for  y^  >  1  and  0  elsewhere,  2Lnd  we  assume  that  along  the  boundary,  for  leu'ge  s 

(6.3.3)  log  s  =  -I-  2a_i  log  -I-  2  ^oo  -I-  oi  a~^  -1-  -t- . .  .^ 

Then 

(6.3.4)  Vo  =  =  ^-‘(a^)"-'  "  , 

v27r 


(6.3.5)  uo  «  s~*(logs)“-‘e“'>v^  , 
and 

(6.3.6)  Voy  =  -S-^a(i>{a)  =  _s-3/2(c^2)«-i  +  l/2_^ga,a-’-Ha,c»-‘  +  ... 

v27r 
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(6.3.7) 


In  the  meantime 


Kiy)  =  -^iV  ^(21ogy2  -  2) 

=  ny~'‘(6logy^  -  lo) 
M'\y)  =  -»'iy~M24logy2-52) 
^i^^(y)  =  riy~®(l201ogy2  -  308) 
h[^\y)  =  -r,y-‘'(7201ogy2  -  2088) 


and  so 

vi  =  ri^ 

(6.3.8) 

and 

Vlj,  =  -riS~^/^Q;“ 

(6.3.9) 

We  would  like  to  have, 


*  *o  ^  log(sa'^)  +  ■* 

1  - 


61og(sa*)  —  lOj 

‘'l201og(ja“)-308j  +  ...| 


21og(AO;^)  —  2|  +  |24log(sa^)  —  52 


1 

2‘ 

1 

n“ 


-6 


720log(sQ2)  -  2088 


+  . 


(6.3.10)  roUo  +  vi^cis  ^ 
and 

(6.3.11)  roVoy  +  VlyWO. 

Since  wi,  «  — 2ris“®/^(logs)”^/^,  (6.3.11)  requires  a_i  +  1/2  =  —1/2  or 

(6.3.12)  a_i  =  -1 
and 

(6.3.13)  roc“®/v^  =  — 2ri 
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But  then  (6.3.10)  requires 


(6.3.14) 
Note  that 


ri  =  c, 


los;(sa^)  =  —  logo^  +  2ao  +  2aiO  ^  +  2o20  ^  +  . . . 

and  thxis,  applying  (4.55)  and  (4.56) 

..  =  /.,{.>'=  a)  +  +... 

vi  =  Cis~^  {  ^  ~  2aia~*  +  2a20 

+  [30'^  +  (—3  log  +  6ao  —  5)o~'*  +  6ai  o"® 

ll5a-‘‘  +  (-15 log +  30ao  -  77/2)a-^ 


(6.3.15) 

and 


+ 


+  1050-”  + 


ui»  =  Kis^^^a)  +  + 


viy  =  —Cis~^^^a~ 


(6.3.16) 
Now 

(6.3.17) 
and 

(6.3.18) 


2  +  (-21ogo^  +4ao  -  2)a“^  +4aiQ"'*  +4020 
12o-2  +  (-121ogQ2  +2400  -26)o-'‘  +24oio-®J 
OOo-'*  +  (-901ogo2  ^  _  261)o-®j  +  840o-®  +  . . .  | 

rouo  =  — 2cis-*|o-^  +  oio”^  +  ^02  +  ■*"  "*■  +  •  •  •  | 

roWoy  ~  2CiS  ^^^O  '/l+OiO  ^+^02  +  ®  ^ 


( 


+  I  03  +  ^102  + 


rovoy  +  Vly  Cl5 


b'i 


2ai  +  2  log  o*  —  4ao  +2  —  12 


+ 


...} 


rovo  +  vi— cis  *  Rs  cis  ^  |o-^  |— 2  —  logo^  +  2ao  +  3  +...|. 
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A 


We  let 


(6.3.19)  h2iy)=y  ^|r2i  logflogy^]  +  r22| 

which  gives  rise  to  t;2  =  ^2  +  sh^2^/2  +  . . .  and  V2y  =  h'^  +  f2  +  . . .,  or 


V2 


W  S 


r2i  log 


o''  —  logo"*  +  2ao  +  2aio~^  +  . . . 


+  r22  + 


] 


—  s  *a  ^|r2i  logo^  +  r22  +  »'2ia  ^  ^-logo^  +  2%  +  2aiO  ^  +  . .  .^  + 

(6.3.20)  V2  «  logo^  +  r22  +  r2io"^^-logo^  +  2ao  +  . .  .^  +  . . .  | 


and 

(6.3.21)  V2y 

from  which  it  follows  that 


=  —s  ^/^o  *|2r2i  logo^  +  2r22  +  •  •  I 


2ci  —  2r2i  =  0 
Ci(2ai  —  4ao  —  10)  —  2r22  =  0 
-Cl  +  r2i  =  0 

Ci(— 2  +  2ao  +  3)  +  r22  =  0 


or 

(6.3.22) 

r2i  =  Cl 

(6.3.23) 

r22  =  -Ci(2ao  +  1)  =  Ci(aj  -  2ao  -  5) 

and 

(6.3.24) 

Oi  =  4  . 

Thus  our  current  approximations  involve  ein  unknown  parameter  oq  .  We  have 


(6.3.25) 


log  s  =  o^  —  2  log  o^  +  2ao  +  So  ^  +  . . . 
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and 


U  = 


2\/^c  ““^(a) 

(6.3.26)  +  y  ^logz^  +  log(logz^)  -  (2ao  +  1)  d2  +  ...| 


The  integral  in  the  expression  above  does  not  converge  because  for  small  z  ,  log(log  ) 
is  not  defined.  For  asymptotic  purposes  the  value  of  the  integrand  in  a  boimded  range 
of  z  values  is  unimportant  and  the  result  would  be  more  meaningful  if  the  integrand 
were  replaced  by  zero  over  the  range  \z\  <  1.  ChemofF  (1965a)  has  a  detailed  discussion, 
for  the  problem  of  deciding  the  sign  of  the  mean,  of  why  we  should  expect  asymptotic 
approximations  of  the  form 


(6.3.27)  u  =  s^^^|ro^(a)  +  TqS  ^^^Q^(a)  +  J  <i>{- ^ )h{z)dz  +  . .  .^ 

and  why  they  lead  to  asymptotic  approximations.  In  the  symmetric  case  r*  =  0. 


While  we  still  lack  control  of  uq,  which  is  crucial  in  u,  it  plays  a  less  prominent  role 
in  the  expansion  for  the  optimal  boundary  which  can  be  written 


(6.3.28)  =  log  s  +  2  log(log s)  -  2ao  +  (log  s)“*  [2  log(log  s)  -  8  -  2oo]  +  . . . 

and  for  which  the  prominent  terms  do  not  depend  on  s.  From  one  point  of  view  this 
states  that  for  large  s  one  should  stop  and  reject  bioequivalence  when  the  P  value 
2^-y/y/s)  is  less  than  the  nominal  significance  level. 


2$(-d) 


2<i(d) 


a 


Since  Co  is  not  known,  this  nominal  significance  level  is  as  yet  unspecified,  but  the  order 
of  magnitude  is  clesir. 


Finally,  one  word  of  caution.  At  this  point  it  is  not  cleEir  that  the  asymptotic  expan¬ 
sion  derived  by  considering  more  corrections  will  be  as  simple  as  (6.3.3).  It  may  require 
replacing  some  of  the  coefficients  01,02,...  by  terms  involving  logo^  or  even  log(loga^). 


6.4  Large  s,  Problem  2. 

There  is  no  need  to  elaborate  here.  This  problem  reduces  to  d2y  =  C2s“*  -4-  1  for 
s  >  1,  and  hence  has  the  seime  expansion  as  in  Section  6.3.  The  only  difference  is  that  ci 
is  replaced  by  C2  and  we  have  a  different  unknown  value  for  oq. 
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